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properties of one of these models, first introduced in 1231 . with more details. By 
introducing a canonical ensemble we calculate the canonical partition function of 
this model exactly. Using the Yang-Lee theory of phase transitions we spot a 
second-order phase transition from a power-law to a jammed phase. The density 
profile of particles in each phase has also been studied. A simple generalization of 
this model in which both the left and the right boundaries are open has also been 
introduced. It is shown that double shock structures may evolve in the system 
under certain conditions. 
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1. Introduction 

The occurrence of non-equilibrium phase transitions is one of the most interesting 
properties of one-dimensional out-of-equilibrium models EI- In open boundary 
systems, changing the boundary parameters may result in special kinds of non- 
equilibrium phase transitions which arc called the boundary induced phase transitions. 
The existence of shocks is another interesting property of these lattice models. A 
shock is defined as a discontinuity between a hight-density and a low-density region 
in the density profile of particles on the lattice. The shocks have been studied both at 
the macroscopic level using the hydrodynamic equations and also at the microscopic 
level by solving the master equation of spatially discrete models. For example, at 
the macroscopic scale the shocks may appear in the solutions of the Burgers equation 
with particles conversation 0. At the microscopic level the Burgers equation can be 
described by the Asymmetric Simple Exclusion Process (ASEP) .IJ. The ASEP cotains 
one class of particles which can be injected and extracted from the boundaries of a 
one-dimensional chain while hopping in the bulk with asymmetric rates. This model 
has been studied widely during the last decade. It is known that by changing the 
injection and extraction rates, the ASEP undergoes several boundary induced phase 
transitions . The shocks may also appear in the density profile of the particles with 
the ASEP dynamics for certain values of the reaction rates. The same phenomenon 
might happen when the fixed number of such particles hop on a ring in the presence 
of a slow particle called an impurity 

Recently, it has been shown that among one-dimensional reaction-diffusion systems 
with open boundaries there are three families of models with a common property: The 
time evolution equation of a Bernoulli shock measure developed by the Hamiltonian of 
each of these models is quite similar to the time evolution equation of a single-particle 
random walk on a lattice with reflecting boundaries |22|. This actually takes place 
provided that some constraints on the reaction rates of these systems are fulfilled. 
These models are the Asymmetric Simple Exclusion Process (ASEP), the Branching- 
Coalescing Random Walk (BCRW) and the Asymmetric Kawasaki-Glauber process 
(AKGP). It is known that in an infinite system, and under some constraints, single 
shock or even consecutive multiple shocks may evolve in the ASEP JOJ-^S]. In the 
latter case the time evolution equation for the multiple shock measure generated by 
the Hamiltonian of the system is similar to the time evolution equation of n random 
walkers. In the steady state the ASEP's probability distribution function can be 
written in terms of superposition of shocks. From the Matrix Product Formalism 
(MPF) point of view, in which the steady state weights are written in terms of the 
expectation values of product of non-commuting operators 5 , n consecutive shocks 
correspond to an (n + l)-dimensional representation of the ASEP's quadratic algebra 

[El El 

The question that might arise is that whether or not the steady state of the BCRW and 
that of the AKGP can also be represented by finite dimensional representation of their 
quadratic algebras since their steady states can be written in terms of interactions of 
single shock. We will deal with this question in the first part of this paper. We will 
show that as far as a single shock structure is concerned the steady state of the BCRW 
and that of the AKGP can be obtained from two-dimensional representations of their 
quadratic algebras. 

In |2.'j| we have studied a branching-coalescing model with open boundary in which the 
particles diffuse, coagulate and decoagulate on a one-dimensional chain. The particles 
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are also allowed to enter and leave the chain only from the left boundary. This model 
had already been studied with reflecting boundaries in ^JEH- We have shown that 
the steady state of this open boundary model can be obtained using the MPF provided 
that the reaction rates satisfy a constraint. In this case its quadratic algebra has a 
two-dimensional representation, while in the case of reflecting boundary conditions 
the quadratic algebra of the model has a four-dimensional representation. Recent 
investigations show that some of the classical concepts and theories in the equilibrium 
statistical mechanics, such as the Yang-Lee theory of equilibrium phase transitions 
|2(J|. can be successfully applied to the out-of-equilibrium systems. There are several 
examples from the driven diffusive systems to the directed percolation models which 
confirm the applicability of this theory to the non-equilibrium systems HS]-j3n|. We 
have also studied the phase transition of our model using the Yang-Lee theory by 
studying the roots of its grand canonical partition function |23| . These zeros lie on a 
circle and intersect the positive real axis at an angle -|. According to the Yang-Lee 
theory this reveals a first-order phase transition As we will see this model is an 
special case of the BCRW; therefore, we expect to see the shock profiles in its steady 
state. In order to see the shock profiles we will confine the number of particles by 
working in a canonical ensemble. The sum over the steady state weights with fixed 
number of particles is defined as the canonical partition function of this model. The 
canonical partition function can be calculated exactly using the MPF. It turns out 
that the zeros of this function in the complex plan of one of the parameters lie on a 
curve which intersects the positive real axis at an angle ?. At the same time their 
density vanishes near the critical point. It indicates a second-order phase transition. 
In the canonical ensemble one can also calculate the density profile of particles on 
the chain again using the MPF. Our calculations show that in one phase the density 
profile of particle has an exponential behavior while in the other phase, as we expect, 
it is a shock. 

Another question that might arise is that under what conditions multiple shocks may 
evolve in this system. We will show that in our model if we let the particles to enter 
and leave the chain from the both ends of the chain, one can see the double shock 
structures given that the injection and extraction rates satisfy some constraints. In 
this case the time evolution equation of the shock measure is similar to that of two 
random walkers. This phenomenon has already been observed in ASEP \22\ I31ll3*2| . 
This paper is organized as follows: In the second section we will briefly review the 
basic concepts of the MPF. In the third section we will show that the quadratic algebra 
of the BCRW, as a generalized form of our branching-coalescing model, has a two- 
dimensional representation provided that the same constraints introduced in |22| are 
fulfilled. In this section the two-dimensional representation of the AKGP will also 
be introduced. In the forth section we will switch to our model proposed in and 
find its canonical partition function using the MPF. The numerical estimates of the 
canonical partition function zeros, as a function of a complex reaction rate, will be 
obtained. The line of zeros and also their density will be calculated exactly. Finally, we 
will calculate the density profile of particles in each phase. The fifth section is devoted 
to the study of double shock structures in a general branching-coalescing model with 
open boundaries. In the last section we will summarize our results. 
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2. The matrix product formalism 

In this section we review the basic ideas of the Matrix Product Formalism (MPF) 
first introduced in 0. We will consider only one-dimensional stochastic systems 
with nearest neighbors interactions and open boundaries. Systems with long range 
interactions and open boundaries have also been studied in related literatures |33| . 
Let us define P(C\t) as the probability distribution of any configuration C of a 
Markovian interacting particle system at the time t. The time evolution of P(C; t) 
can be written as a Schrodinger equation in imaginary time 

j t P{C;t)=HP{C;t) (1) 

in which H is a stochastic Hamiltonian whose matrix elements are the transition rates 
between different configurations. For the one-dimensional systems of length L which 
consist of one species of particles the Hamiltonian H has the following general form 

L-l 

H=^2 h k,k+i +h 1 + h L (2) 
fc=i 

in which 



ifc,fc+i 



: 1®^ ®h® l®( L - k -V , hi = h® <g> 2®(i-i) , h L = <g) 



where X is a 2 x 2 identity matrix, /i is a 4 x 4 matrix for the bulk interactions and h"> 
(h^) is a 2 x 2 matrix for particle input and output from the left (right) boundary. 
In a basis (00, 01, 10, 11) the bulk Hamiltonian along with the boundary Hamiltonians 
are given by 



/ . 


Wi2 


W13 


W U \ 


W21 




W23 




W31 


W 3 2 




W34 


\ mi 


W i2 


W43 


, ■ J 



(3) 

"<"=(": a ).» w - (1 -i 

Requiring the conservation of probability, the diagonal terms of matrices should be 
the negative of the sum of transition rates in each column. In the most general form 
the interaction rates are 

Diffusion to the left and right W32, W23 

Coalescence to the left and right W34, W24 

Branching to the left and right u> 43 , W42 

Death to the left and right Wi3 ; W12 , 

Birth to the left and right W31, W21 

Pair Annihilation and Creation wu, W41 

Injection and Extraction at the first site a, 7 

Injection and Extraction at the last site S, (3. 

In the stationary state we have HP*(C) = 0. According to the MPF the stationary 
probability distribution P*(C) is assumed to be of the form 

1 L 

P*(C) = P*(ti,---,t l ) = —{W\ Y[(nD + (1 - n)E)\V) (5) 
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where Tj = if the site i is empty and t,; = 1 if it is occupied by a particle. The function 
in JSJ is a normalization factor and can be obtained easily using the normalization 
condition ^ c P*(C) = 1. The operators D and E stand for the presence of a particle 
and a vacancy at each site. These operators besides the vectors (W\ and |V) should 
satisfy the following algebra 



E 
D 



(W\ 



E 
D 



E 
D 



h (r 



E 
D 



\V) 



-(W\ 

E 
D 



E 
D 

E 
D 

\V). 



E 
D 



E 
D 



E 
D 



(6) 



The operators E and D are auxiliary operators and do not enter in the calculation of 
physical quantities. Using © and © the quadratic algebra associated with the most 
general reaction-diffusion model can be obtained. By defining C := D + E one can 
easily see that the mean values of the physical quantities can be written in terms of 
the non-commuting operators C and E and the vectors |V) and (W\. For example 
the mean density of particles at site i has the following form: 



(Pi) = 



E,c T i p *( e ) _ (w\c i -\c - E)C L - l \V) 



Similarly, any n-point density correlation function can be written as 

Y,cTnn 2 ---T in P*(C) 



(7) 



(PhPii ■ --Pi n ) 



P*(C) 



(C - E)C L -^\V) 



(C L ) 



(8) 



In order to calculate Q, ©, or any other quantity, there are two possibilities: one 
might find a finite or an infinite-dimensional representation for the quadratic algebra 
of the system or equivalently, one may use the commutation relation between the 
operators to calculate these quantities rigorously . In the next section we will show 
that finite-dimensional representation exists for three families of reaction-diffusion 
systems, given that some constraints are satisfied. 



3. Representation of the quadratic algebra 

In our recent paper |22j exact travelling wave solutions are obtained for three families 
of one-dimensional reaction-diffusion models with open boundaries. These models 
are the Asymmetric Simple Exclusion Process (ASEP), the Branching- Coalescing 
Random Walk (BCRW) and the Asymmetric Kawasaki-Glauber Process (AKGP). 
It has been shown that for these models the stationary measure can be written as a 
linear combination of Bernoulli shock measures defined on a lattice of length L as 

H = ( j ® ( ~ ^ J , < m < L (9) 

provided that some constraints on the reaction rates are fulfilled. In (jSJ) p\ and pi 
are the densities of particles at the left and the right domains of the shock position 
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m respectively. The time evolution of 10 generated by the Hamiltonian of the above- 
mentioned models is given by 

- H\m) = di\m - 1) +d 2 \m + 1) - (di +d 2 )|m) , 0<m<L (10) 

which is a simple single-particle random walk equation for the position of the shock 
m with hopping rate di to the left and d 2 to the right [22j- In this section we 
show that under the same constraints the stationary states of these models can be 
studied using the MPF and their associated quadratic algebras have two-dimensional 
representations. 

For the ASEP the non- vanishing rates in Q are W32, W23, a, /3, 7 and 5. It is known 
that the quadratic algebra of the ASEP has an n-dimensional representation (n > 2), 
provided that the parameters of the model satisfy some constraints 117| . The 
(n + l)-dimensional representations of the quadratic algebra describe the stationary 
linear combination of shock measures with n consecutive shocks. The ASEP has 
already been studied widely in literatures; therefore, will not be considered here (see 
and references therein). 

For the BCRW the non- vanishing parameters in are W34, W24, w± 2 , 1043, W32, 
W23, at, (3 and 7. As we will see, the steady state of this model can be described by 
two-dimensional representation of its quadratic algebra given that 

p UJ42+'!i>43 p (11) 

7 = ^r a + (1 _ p)w?,2 - ^7^43 + pw 34 

in which p\ = p is the density of particles on the left-hand side of the shock. The 
density of particles on the right-hand side of the shock p 2 is zero. 
For the AKGP, the non- vanishing parameters in Q are W12, W13, W42, W43, W32, a 
and /3; however, there is no specific constraint on these parameters in order to have 
the mentioned property. If the particles are allowed only to enter the system from the 
first site with the rate a and leave it from the last site of the chain with the rate (3, 
then we will have p\ — 1 and p 2 = 0. 

We have found that the quadratic algebras of the BCRW and the AKGP have two- 
dimensional representations given that the constraints for the BCRW are fulfilled. 
For en ^ C22 we find 

c =( ^ ° ) > E =( G n Gl2 ) 
c 22 J \ e 22 J 

(12) 



c = 1 c C12 E _ I en ei2 



(13) 



and for c = en = C22 we find 

c 1 

5 J ' ~ ~ V e 22 

The matrix elements are functions of the non-vanishing parameters of each model. 
The explicate form of the matrices C, E, C and E and the vectors \V) and (W\ for 
each model is given in Appendix. Using the representation l|12|) for cu ^ c 2 2 and the 
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definition of the density profile of particles (JJJ one finds in the thermodynamic limit 
(L -> oo) 

en v 2 ei2 kzk 
(1 ) e^, en > c 22 

I \ ) Cn ViCu 

<*>=< n e 22 , «Jie ia -i ( 14 ) 

(1 j e « , Cn < c 22 - 

C22 w 2 c lx 

in which we have defined the characteristic length £ = | ln^-| . Also using (|13|) and 
(J7J) for c = cn = C22 

(Pi) = (1 - ^) + - %r , < x < 1 (15) 
c c c 

where x = -f. As can be seen the density profile of particles has an exponential 
behavior for cn / C22 while it is linear on the coexistence line cn = 022- For Cn > C22 
the density of particles near the left boundary is p± = 1 — ^ while for cn < C22 the 
density of particles near the right boundary is p 2 — 1 — f^- ■ 

Having the two-dimensional representations i|12|) and l|13l) any higher order correlation 
functions can easily be calculated using ||SJ|. For instance, in the thermodynamic limit 
the two-point correlation function is obtained to be 

dn , 2 v 2 dndi2 m , v 2 d 12 d 2 2 
( ) H 5 — e « -I e « ,cii>c 2 2 

m = { T v f\ V1 TT 



d 22 s 2 , w\di\di2 rl , w 1 di 2 d 2 2 _ 

( J H 2 — e 5 "! e 5 i c H < c 22- 

C22 W2Cn W2C11C22 



for cn 7^ C22 and 



/ \ / d iiN2 , d\i{d 2 i - dn) j d 2 2(d 2 2 - dn) i 

( PiPj ) = (— ) + =5 ( z ) + j 2 ( z ) (17) 

for c = cn = C22- In ((TCI) d^'s are defined as Aj = (C — i?)y = rfij and in ((171) d^'s 
are the same elements for the case Cn = 022- 



4. The Yang-Lee zeros of the canonical partition function 

In this section we will focus on some of the steady state properties of our branching- 
coalescing model first studied in We will first calculate the canonical partition 
function of our model with the following non-zero rates: 

W32 = q,W 2 3 = ^,W 34 = <7,^43 = Ag,W 2 4 = \ , W42 = — , a, 7. (18) 

It can easily be verified that this model is an special case of the BCRW. In this model 
we have (3 = and the constraints (|llfl also give p — and a = (~ — q + ~/)A. 
By introducing the canonical ensemble we restrict ourself to the case where the total 
number of particles on the chain M is fixed. As we mentioned, it has been shown that 
the classical Yang-Lee theory of the equilibrium phase transitions can successfully 
be applied to the driven diffusive systems to predict their non-equilibrium phase 
transitions. By finding the roots of the canonical partition function of our model, 
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we will investigate its phase transitions. In order to calculate the canonical partition 
function we will use the quadratic algebra of the model 

[C, C] = [E, E] = 



EC ~ EC ={q + qA + q- x )EC - q(l + A)E 2 - q^C 2 

CE-CE= {q- 1 + q^A + q)CE - q~ l (l + A)E 2 - qC 2 

(W\((a + i)E + E- 7C) = (W\C = , E\V) = C\V) = 
and its representation for q 2 ^ 1 + A 

1 + A \ „ ( 1 A \ ^ ( ^— iA 
q* ) > E =[0 q 2 ( ^ 



(19) 



(20) 



For convenience we will use our new notation in which 7 (instead of (3 in |23p is the 
extraction rate of particles from the left boundary. It can easily be seen from (I2UII 
that the operators C and E have the following properties 

E(C-E) = {C-E) , (C-E) i = A i - 1 (C-E) (21) 

We define the canonical partition function Zl,m as the sum over the steady states 
weights with fixed number of particles M. Using p9 (l -l|21 [l and the approach used in 
IBHH! we obtain 

z "' - 7( i + a")- ,a § <?' c ^-> < 22 » 

where Cf = .,, t! .,, is the binomial coefficient. Let us now examine the zeros of the 
canonical partition function l|22|) as a function of q for fixed L, M, a, 7 and A. Apart 
from the zeros generated by the factor behind the sum in i|22|) , the numerical estimates 
of the roots in the complex-g plane are plotted in Figure (1) for L — 500 and M — 300 
(grey dots). As can be seen from Figure (1) the roots accumulate towards a critical 
point on the positive real-q axis which will be calculated shortly. In order to calculate 
the line of zeros analytically we need the thermodynamic behavior of the canonical 
partition function l|22ll . This can be done using the steepest decent method. We find 

q 2L A M l „ 
P < 1 - '/ 



ZLM { P CfA™ 2 (23) 

7-7— It TTj p>l-q- 2 

in which p = ^r. The critical point is obviously q c — ^l_ p - Defining the extensive 
part of the free energy as 

g= lim y\nZ LM (24) 

L,M— *-OC J_j 

one can calculate the line of zeros from 

Re gx = Re g 2 (25) 
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Figure 1. Numerical estimates (grey dots) for the roots of y"]^_ n M q 



besides the line of zeros (dark line) given by 1261 for L - 
complex-q plane. The exact critical point is q c = J— 



: 500 and M 
1.581. 



A/-1 
L-i-1 
300 in the 



in which g\ and gi are the free energy functions on the left and the right-hand side of 
the critical point [21] • After some calculations using l|23 |l -(|2ft |l we obtain the following 
equation for the line of zeros in the thermodynamic limit: 



x 2 + y 2 



1) 



{2xyf 



P" 



(26) 



[{x 2 - y 2 

in which we have defined x := Re(q) and y := Im{q). This function is plotted in 
Figure (1) for L = 500, M = 300 and p = ™ = 0.6 (dark line). It is seen that the line 
of zeros lays on the numerical estimates of the zeros. The function (|26|l intersects the 
positive real-q axis at the critical point x c — yj= • We can also predict the order of 
transition by fining the density of roots near the critical point. For < y <C 1 and 



(x- 



= ) <C 1 wc find that the line of zeros is actually a straight line y 



with the slope j. The density of zeros on this line can be computed from 

^s) = ^\-l~Im( gi -g 2 )\ (27) 
2n os 

where s is the distance from the transition point It turns out that the density 

of zeros as a function of the distance in the vertical direction is proportional to y so 
in y — > limit the density of zeros vanishes. This expresses a second-order phase 
transition at the critical point. 

In what follows we will calculate the density profile of particles on the chain in each 
phase using the canonical ensemble. In the canonical ensemble, the density of particles 
at site i can be written as 



{Pi 



(28) 



Note that the sums are over the configurations with fixed number of particles M. By 
using 1)19 (1 - l|2i p and after some calculations we obtain 



L-fc-2 



+ a 



M-l 2(L- 
2 U 



)0(M < i) 



2^k=0 1 ^L-k- 



(29) 



in which 0(- ■ ■) is the ordinary Heaviside function. In Figure (2) we have plotted (|29|l 
for L = 500, M — 300 and two values of q. With this choice of L and M the critical 
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Figure 2. Plot of the equation E^J for L = 500, M = 300 and two values of q 
as it is shown. Below the transition point (q c = 1.581) the density profile has an 
exponential behavior while it is an error function above the critical point. 



point is q c 



1.581. As can be seen for q = 1.5 < q c the density profile of 



particles is constant through the bulk of the chain equal to p = = 0.6 and drops to 
zeros only near the right boundary. For q = 2.5 > q c the density profile of particles is 
a shock. The high-density region (puigh = 1 — q 2 = 0.84) extending from i = 1 to 



_ p_ 



sL is separated by a sharp interface from a low-density region (plo 



0). 



Using the steepest decent method one can easily show that in the power-law phase 
(q < q c ) the density profile of particles (pj) has the following analytical form 



L-i+l . 
5 



(Pi) = P(l ■ 

in which £ is the characteristic correlation length equal to £ 



(30) 



ln(? 2 (l-p))|- 



5. Generalization: Double shock structures 

In our non-conserving driven diffusive model 23 only the left boundary was assumed 
to be open so that the particles could enter and leave the system from the first site of 
the chain. In this case the time evolution of the single shock measure was equivalent 
to that of a random walker on a one-dimensional lattice with reflecting boundaries. 
In this section we will introduce a simple generalization of this model in which the 
particles are allowed to enter and leave the system from both ends of the chain. Our 
main objective will be to investigate the possibility of the existence of double shock 
structures and the way that they evolve in time. Defining a double shock measure, 
we will see that on special manifold of parameters its time evolution is equivalent to 
that of two random walkers on a one-dimensional lattice with reflecting boundaries. 
Let us assume that apart from the interactions given by (|18|) . the following reactions 
take place at the right boundary 

(3D 



injection at the last site with rate 
extraction at the last site with rate f3. 



On a chain of length L we define an uncorrelated double shock measure as 

/ , \ ®m / \ Xri m 1 / \ ®L — n+1 

\ mt n)=( 1 -») J 1 '/ 2 ) 4 ;M (32) 
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Pi 1 | 

P2 I 1 

1 m n L 

Figure 3. Density profile of a double shock measure 



in which m + 1 < n, < m < L and 1 < n < L + 1. We have sketched the density- 
profile of a sample double shock measure in Figure (3) . It can easily be verified that the 
time evolution of the double shock measure is equivalent to that of two random 
walkers provided that 

Pl = jf^, P2 = 0, a = (q- 1 -q + j)A, S = (q - q' 1 + 0) A.(33) 

Under these conditions the time evolution equation of l|32(l generated by the 
Hamiltonian of our generalized model will be given by the following equations 

—H\m, n) = q\m — 1, n) + q (1 + A)|m + 1, n) + q^ 1 \m, n + 1} + 
q(l + A)\m,n- 1) - (2 + A)^ 1 +q)\m,n), 
for m =/= 0,n =/= L + l,n m + 1 

-H\0,n) = a(^)\l,n) +q- 1 \Q,n+l) + q(l + A)|0,n- 1)- 



(a(^)+q-'tq(l + A))\0, 
for n^l,£ + l 



n > 



-H\m, L+l)= + A)\m + 1, L + 1) + q\m - 1, L + 1>+ 

5(±±A)\m, L) - (5(±±A) + q + q -^l + A))\m, L + 1), 
for m/0,L 

-H\0,L+l)=a(^)\l,L + l)+S(^)\0,L)- 
(°H^l+A) | 0)£+1) 

— H\m,n} = for n = m + l. 

The last equation indicates that the model has a trivial steady states in which the 
density of particles at all sites is equal to pi = jA^ . This is associated with a one- 
dimensional representation of the quadratic algebra of this model. As can be seen from 
<|34[1 the shock fronts move like two biased random walkers with different hopping rates. 
The shocks are also reflected from the boundaries. In the most general form it can 
easily be seen that for an open boundary model with the following non- vanishing rates 

W23, W 32 , W 2 4, W42, ^34, W43, a, (3, 7, 8 (35) 



(34) 
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the above-mentioned approach can be used to show that the double shock structures 
of kind i|32|) exist if the reaction rates satisfy the following constraints: 

p _ W42+W43 _ Q 



_ «'4 2 (u'24 + t"34) W i 3(w 2 4.+W 3 4,) 

32 _ W i2 +W i3 ' 23 _ IO42+MI43 

n _ W42+W43 f _ 2W4 2 (W 2 4+W34,) , (m 2 4 + tU34)(tt)34+m 4 2) . \ 

V 2 ^ WAO+WA3 WOA +W3A +WAO+WA3 ' ' 



(36) 



W24+W3S, \ ^ W42+W43 U) 2 4+U) 3 4+U)42+J«43 

J W42+W43 / , 2M)4 2 (tU24+tU34) _ (wgj + W34 ) (t»34 +W42 ) . g\ 

MI24+W34 ^ 24 W42+W43 M>24+M>34+M>42+1A'43 

In this case the time evolution equations for |m, n) given by l|32|) will be quite similar 
to (j34|) . One can readily check that for our generalized model, defined by l|18|l and 
13 It . the conditions l|36(l reduce to l|33|l . The steady state of the system can be 
obtained from the superposition of double shock measures |34j : however, the system 
has a trivial steady state in which the particles occupy the sites of the chain with 
the uniform probability p\ and as before this is associated with a one-dimensional 
representation of its quadratic algebra. 



6. Concluding remarks 

Our calculations in this paper can be divided into three parts. In the first part of 
the paper we have shown that using the MPF the steady state of the ASEP, BCRW 
and AKGP can be expressed by two-dimensional representations of their quadratic 
algebras provided that some constraints hold on the reaction rates. These constraints 
introduce specific manifold of parameters. On these manifolds the stationary measures 
are given by superpositions of Bernoulli shock measures. Having the explicit form of 
the two-dimensional representations several physical quantities, such as correlation 
functions, can easily be calculated. 

The BCRW is a generalized form of the model that we had studied in The 
second part of our calculations is devoted to the detailed study of the phase transition 
and also the shock formation in this model as an special case of the BCRW. We have 
introduced a canonical ensemble in which the number of particles on the chain is fixed. 
Using the properties of the quadratic algebra the canonical partition function of the 
model is calculated exactly. The thermodynamic behavior of this function reveals a 
phase transition which depends on q and p i.e. one of the reaction rates and the 
density of particles on the chain. In order to study this phase transition we have used 
the classical Yang-Lee theory. Our calculations show that the Yang-Lee zeros of the 
canonical partition function as a function of the complex variable q lie on a curve 
which intersect the positive real-q axis at an angle j at the critical point q c = ■ 
According to the Yang-Lee theory of phase transitions, this is reminiscent of a second- 
order phase transition. The nature of each phase can be studied by calculating the 
density profile of the particles in the steady state. For q < q c the density of particles 
is constant everywhere on the chain except near the right boundary where it drops to 
zero exponentially with the correlation length £ = \q 2 (l — while it has a shock 

structure for q > q c . 

In the third part of our calculations we have shown that double shock structures can 
evolve in a general branching-coalescing model with open boundaries. In this case the 
boundary parameters lie on specific manifolds of parameters determined by (|36|) . A 
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detailed study of the double shock properties will be given elsewhere |34j . 
The exact solution of the BCRW without any constraints on its reaction rates still 
remains an open question. The exact steady state of the model defined by the non- 
vanishing rates (|35|) under the conditions (|36|l from the MPF point of view is under 
consideration |34| . 



Appendix 

It turns out that the representation of the quadratic algebra of the BCRW in the case 
en 7^ C22 (see JEJl) is 

C=(& ° ) E=( W ™ A ^ 

(1 - p)w 32 + pw 3i J ' \ (1 - p)w 32 + pw u J' 

E j -^43(W23 ~ W 32 (l -p)~ W 3iP ) "Af23i ^ ^ ^ ^ 

(A(l+-2-) \ 
i3+w 32 (l- P )lw 3 4P-io 23 ),lW\ = ( y, T ^Ti^lT^l !T,lY 
J ] \ A(q + UJ 2 3 ^ T +p(l«32(l-p)+-U'34P)) ' / 

and in the case en = C22 (see 1(15)) it is 

^ _ / w 32 (l - p) + w 3 4P A 

V w 32 (l - p) + w 34 p 

E = ( i 1 - P)Ml - P) + w 34 p) p{w32{1 Z V p)+Waip) 

\ ^32(1 - Pi) + W34P J ' 

E = ( ^ (W32(1 7> + W34P)2 J ) , C = , |F) = ( * ) , (W\ = («*, 1 
in which 

_ V0+W32(:n+^l3)p-(w 3 2-W34)(ri+\p)p 2 

1 ~ p 2 (lU 32 (l-p)+U'34p) 2 (/3+p(lU32(l-p)+-U)34P)) 

w _ ap 2 (-U->32(l-p)+^34P) 2 

1 JJQ + Ap(lU32(l-p)+lil34p)(a(l-p)+p 2 (M)32(l-p)+U'34p)) 

and A and 77 are arbitrary parameters. The representation of the quadratic algebra 
for the AKGP is also found to be 



W 13 J ' V W 13 / V 

X{w i3 -w 13 +l3) 



and 



I / I V A(wi3-UJ43+Q) ' 



IV) = ( f ),W = (j^.i)- 

for the case en / C22 and Cn = C22 respectively. 
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